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We study the two-dimensional Bose-Hubbard model in the presence of a 3-body interaction term, 
both at a mean field level and via quantum Monte Carlo simulations. The 3-body term is tuned 
by coupling the triply occupied states to a trapped universal trimer. We find that, for sufficiently 
attractive 3-body interaction the n = 2 Mott lobe disappears and the system displays first order 
phase transitions separating the n = 1 from the n = 3 lobes, and the n = 1 and n = 3 Mott 
insulator from the superfluid. We have also analyzed the effect of finite temperature and found that 
transitions are still of first order at temperatures T ^ J. 



The Bose Hubbard (BH) model and its continuous 
superfluid (SF)-Mott-insulator (MI) transition represent 
one of the paradigmatic examples of strongly interacting 
many-body physics in lattice structures [l] . The unprece- 
dented control over ultra-cold atoms in optical lattices 
not only allows for a clean experimental realization of the 
BH model [2j|3], but also the exploration of a panoply of 
quantum effects beyond the standard BH model, such as 
disordered induced SF-insulator transitions [4], and the 
reduction of superfluid coherence in mixtures [5]-[7] . From 
the theoretical standpoint, these lattice systems are de- 
scribed by effective low energy extensions to the standard 
BH model which have been shown to give rise to richer 
phase diagrams beyond the SF and MI transition [8HTT] . 
One key element for such impressive progress is the pos- 
sibility of tuning two-body interactions by using Fesh- 
bach resonances or changing the strength of the lattice 
confinement. More recently, effective multi-body inter- 
actions have been experimentally observed [T2', T3^ . The 
question that naturally arises is how these interactions af- 
fect the many-body behavior. Topological phases such as 
fractional quantum Hall states appear as ground states 
to model Hamiltonians with strong three-body interac- 
tions while exotic quantum phases have been predicted 
for bosonic Hamiltonians with many-body interactions, 
such as the ring exchange model [15] . An important 
first step in realizing these models using ultra-cold atoms 
was the recognition that strong three-body losses lead to 
an effective hard core three-body interaction that can be 
used, for instance, to stabilize the BH model with attrac- 
tive two-body interactions ^6j. Under these conditions, 
the system can undergo a first order MI to SF transi- 
tion in the presence of strong pairing interactions [17]. 
Despite these recent studies, lattice systems with three- 
body interactions remain largely unexplored. 

In this Letter, we analyze how the many-body physics 
of the BH model is affected by the presence of local and 
tunable three-body interactions. We propose a mech- 



anism for engineering a 3-body onsite interaction term, 
[/a, which is controlled by an external RF pulse that cou- 
ples the triply occupied state with a three-body bound 
state associated with an excited hyperfine state. This 
local three-body interaction only affects triply occupied 
sites leading to a modified BH Hamiltonian 

H = -JY^ alaj^^l^niirii - 1) ^Sn^^sUs - Itrii] (1) 

where aj(ai) is the bosonic creation (annihilation) oper- 
ator. Hi = a\ai^ J is the hopping matrix element, U is 
the 2-body onsite interaction, U^, is the 3-body onsite in- 
teraction, /i is the chemical potential, and (i, j) denote 
sum over nearest neighbor sites only [see Fig.[TJa)]. Note 
that the three-body interaction considered here is differ- 
ent from the more conventional interaction of the form 
/73n(n — l)(n — 2)/6. Then, we analyze the phase diagram 
of the modified BH Hamiltonian focussing particularly in 
the > 0, [/a < regime. A mean-field Gutzwiller ap- 
proach shows the standard BH phase diagram is modified 
by the presence of attractive interactions. In particular 
for I/73I > [Fig. 2 (d)] the n = 2 (n = 4) lobe disappears 
as double (quadruple) occupancy becomes energetically 
unfavorable. In this regime, we find that, in contrast to 
the standard BH model, a transition between the n = 1 
and n = 3 lobes can occur at finite hopping and that the 
transition between the lobes becomes first order. Such 
transition line extends up to a triple point from which 
two first order transition lines, separating the n = 1 MI- 
SF and n = 3 MI-SF phases, depart. Quantum Monte 
Carlo (QMC) simulations [TF confirm the existence of a 
first order phase transition and provide quantitative pre- 
dictions of the phase diagram in two dimensions and par- 
ticular values of U and Us. Finally, we briefly discussed 
finite temperature effects and experimental signature of 
the first order transition. 

Tunable three-body interaction: To achieve a separately 
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FIG. 1. Schematic representations of (a) the terms in the 
Bose-Hubbard Hamiltonian considered here [Eq. ([T])] and (b) 
the RF field tuned three-body onsite interaction are shown. In 
(b) red circles represent bosons in the lowest hyperfine state, 
6, while blue circles represent atoms in the excited hyperfine 
state, X. 



tunable, on-site, three-body interaction we envision a sys- 
tem in which a universal three-body bound state, at- 
tached to an excited hyperfine threshold, is coupled to 
the identical boson ground state through an external RF 
field, (shown schematically in Fig.jljb)). In this scheme, 
identical bosons in two different hyperfine states (labeled 
b and x for the lowest state and an excited state respec- 
tively) sit on a single site which we model as an isotropic 
oscillator with oscillator length I ho- For a universal three- 
body bound state to form in an excited hyperfine state, 
we assume that < a^ij <C Iho < cibx where a^j is the 
s-wave scattering length between particles in states i and 
j. In this situation, universal three-body Efimov states 
form attached to the bbx three-body hyperfine state [19 . 

The binding energy of such trimer states is determined 
by the three-body parameter, Rq. We envision a system 
where Rq is such that an Efimov state is weakly bound 
with binding energies (in units of h = 1) Ei^ ^ u where 
00 is the effective trapping frequency of the lattice site. 
Under these conditions, the three-body energy for two 
identical bosons with a third distinguishable particle {bbx 
internal configuration) is lower than the two-body energy 
{bx configuration) [20 . Further, in our initial numerical 
studies, we have found that for each Efimov state in this 
scenario, there is a single four-body state {bbbx configu- 
ration) bound below the trimer state. This is in contrast 
to the more commonly considered case of four identi- 
cal bosons in free space in which there are both a deeply 
bound and a weakly bound tetramer associated with each 
Efimov state [21 . The energies for two-, three-, and four- 
body states are shown schematically in Fig. [TJb). 

The above scenario has several benefits. First, the 



large spatial extent of a weakly bound three-body state 
can give good overlap with the ground, identical boson 
state. Second, we can treat a lattice site approximately 
as an isolated two-state systems with energies and Rabi 
frequencies determined by the number of particles on the 
site. Here we focus on two, three, and four-body occu- 
pation, but the general concept carries to higher occupa- 
tion numbers as well [22j . A problem with using Efimov 
states to tune three-body interactions is their generally 
short lifetimes (on the order of 10/is [23]). To avoid this, 
we consider the large detuning limit, ^ ^ O3, where S is 
the detuning of the RF field from the three-body transi- 
tion energy, which can be either positive (red detuned) 
or negative (blue detuned), and is the Rabi frequency 
for the n-body association process. In the large detuning 
limit the effective lifetime of the three-body state is en- 
hanced by a factor of and the shift in the effective 
onsite energy is given by Us ^ 

To ensure that only three-body occupation exhibits 
any significant shift we also require that two- and four- 
body states are not significantly effected by the pres- 
ence of the rf field. To achieve this we require that 
ny{A23^2U^S) < \Us\ and nl/{AsA-3U - 5) < \Us\ 
where Aij is the energy difference between excited states 
with i and j particles. The first of these inequalities 
relies, primarily, on the detuning from the three-body 
transition being much smaller than the energy splitting 
between two- and three-body excited states. 

Achieving large enough energy separations while main- 
taining a reasonable Rabi frequency may require some 
fine tuning. This is because, the three-body bound state 
becomes smaller in spatial extent as it becomes more 
deeply bound leading to a decrease in Qs. The sec- 
ond inequality is somewhat more readily achieved as the 
four-body binding energy is proportional to the three- 
body energy [21] meaning that there is always an en- 
ergy separation. Further, because the four-body state 
is more deeply bound with respect to the three-body 
binding energy, we can expect that its spatial extent is 
smaller leading to smaller overlaps with the identical bo- 
son state and smaller values of O4. We focus here on 
the scenario in which three-body onsite interactions can 
be tuned without significantly effecting other occupation 
energies. However as long as the deviations in three- 
body onsite energies can be tuned to be comparable with 
two-body energies, while four-body energy deviations are 
small, we expect that qualitatively, the behaviors exam- 
ined here will persist. 

Our initial investigations indicate that using the above 
scheme, the three body interaction might be tuned to be 
comparable with the two-body interaction, i.e. \Us\ ~ U. 
We present the above idea here as a plausibility ar- 
gument for a tunable, three-body, onsite interaction. 
A more detailed study to determine the effects of rf 
coupling to an excited three-body state is left for future 
investigations [24 . However, direct rf association of 
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FIG. 2. (Color online) Mean field phase diagram in the fi/U 
vs. J/U plane for \U3\ = (a), \U3\ = U/2 (b), \U3\ = U (c), 
and \U3\ = -3U/2 (d). When < \U3\ < U the n = 2 and 
n = 4 lobes visibly shrink in favor of the n = 3 lobe, until they 
completely disappear at | [/a | = [/ (c) . For | [/s | > U the n = 3 
lobe begins to overlap with the n = 1 and n = 4 lobes and 
a direct phase transition between MI lobes becomes possible. 
The dotted rectangle in (c) highlights the region examined in 
detail in Fig. |3] 



universal trimer states has already been demonstrated 
in ultra-cold, three-component, Fermi gases [25 lending 
credibility to the experimental accessibility of this model. 

Ground state properties of the modified BH model: 
mean field analysis: We will now use Gutzwiller mean 
field theory to study the modified BH Hamiltonian de- 
scribed by Eq. ([T]). 

The Gutzwiller mean field theory is constructed by re- 
placing the full Hamiltonian by an effective local Hamil- 
tonian subject to a self-consistency condition. We intro- 
duce the superfluid order parameter 7/^ = (aj) = (a^) and 
the Gutzwiller wavefunction \G) = H^q(^^q /n*^ |n^)), 
so that the effective Hamiltonian for a translationally in- 
variant system, i.e. fn^ = Z^, takes the form, 

E[ij] = -JzijJ2y^^^ifn+lfn + C-C.)+zJxP^ + E^ (2) 

n 

where z is the coordination number, fn are variational 
parameters, and = ^n{n—l)-\-Sn,3Us—jiin. The prob- 
lem is now reduced to determining the set of coefficients 
{fn} which minimize E[iIj] and satisfy the normalization 
condition {G\G) = Y.n \fn? = 1 [aEHEl!. 

Figure 2 shows the ground state phase diagram of 
model ([T]) at different values of Us. As we increase the 
magnitude of \Us\ from to /7 the n = 2 and n = 4 Mott 
lobes shrink considerably while the n = 3 lobe increases 
in size [as seen in Fig. [2|b)]. In particular, for Us ^ —U 



the n = 2 and n = 4 lobes completely disappear since it 
is now energetically more favorable to have occupation 
number n = 3 [Fig. [2];c)]. This can be easily understood 
in the zero hopping limit. At /i = /ii2 = /7 a doubly 
occupied site has the same energy as a singly occupied 
one. At /i = /ii3 = {3U + Us)/ 2^ instead, a singly 
occupied site has the same energy as a triply occupied 
one. The condition /ii3 < /ii2 sets the Us value for which 
the second lobe disappears, i.e. |/73| > /7 (at |/73| = /7, 
sites with occupation number n = 1,2,3 are degenerate 
in energy for ja = /ii2 = fiis)- Direct transitions from 
MI occupation numbers n = 1 to n = 3 survive at finite 
hopping [Fig. |2|d)] confirmed below using QMC. The 
same argument shows that \Us\ > U also implies the 
disappearance of the 4th lobe. One can easily see that 
upon further increasing Us^ all lobes other than n = 3 
will eventually disappear, e.g. at Us = —3/7, the n = 1 
and n = 5 lobes disappear. 

We have monitored the behaviour of mean field energy 
Eq. Q at fixed /i/U while varying J/U to study the 
order of phase transitions described by model ([T]). The 
formation of double minima structure in the mean field 
energy functional E[iIj] is a signature of first order phase 
transitions. We have observed such double minima 
structure at Us = —1.5/7 for the n = 1 MI-SF and the 
n = 3 MI-SF transitions. The occurrence of first order 
transitions can be understood with a simple argument. 
At fixed small J//7, Us ~ /7, and upon increasing 
(decreasing) fi in order to dope the n = 1 (n = 3) 
MI with particles (holes), double occupancy will be 
suppressed in favor of triply occupied sites. At large 
enough Us such mechanism will eventually prevent a 
gradual addition (subtraction) of particles resulting in 
first rather than second order transitions. Second order 
transitions will be restored at large enough J/U as the 
kinetic energy gain due to hopping of extra particles 
(holes) will again favor a gradual change in density. 

Comparison with QMC simulations: In order to 
confirm the mean- field predictions, we have performed 
QMC simulations on a square lattice of linear size up to 
L = 24 (and L = 30 in certain cases) for selected values 
of J//7, and at /3 = = L/J which corresponds to 
effective zero temperature regime. Figure [3] compares 
the QMC results with the mean-field predictions of the 
phase diagram for Us = —3/7/2. As mentioned above, 
a direct transition from n = 1 MI (in blue) to n = 3 
MI (in pink) survives at finite hopping. This first order 
transition is depicted by the dotted line, while solid 
lines refer to first order transitions from MI to SF. The 
solid and open symbols correspond to QMC predictions 
of the phase boundary with first and continuous phase 
transitions. 

To extract transition points, we have analyzed the 
particle density n as a function of ja. Additionally we 
have performed hysteresis analysis by sweeping back and 



4 




0.04 0.08 0.12 0.16 0.2 

zJ/U 



FIG. 3. First order phase transitions at U3 = —1.5U. Dotted 
line corresponds to the n = 1 MI (blue)-to-n = 3 MI (pink) 
first order mean field transition line. Solid lines refer to mean 
field first order transitions from MI to SF. Solid circles are 
QMC results from hysteretic curves. Dashed lines correspond 
to second order mean field MI-SF transitions. Open squares 
are second order transition points from QMC. Lower (upper) 
inset shows examples of hysteretic behavior for the n = 1 
(n = 3) MI-to-SF transition. 



fourth in chemical potential and calculating the corre- 
sponding particle density. The hysteretic behavior of the 
system along the phase boundaries (solid symbols) fur- 
ther confirms that these are first order transitions. We 
show two examples of such curves for the n = 1 MI-SF 
and n = 3 MI-SF transitions, in the lower and upper 
inset respectively. Based on the energy argument previ- 
ously discussed, we expect the phase transition to become 
continuous as J/U is increased. Indeed, larger kinetic en- 
ergy will favor formation of particle/hole excitations on 
top of the MI. The energy gain due to hopping of the lat- 
ter will compete with the attractive 3-body interaction 
and will eventually restore the second order MI-SF tran- 
sition driven by addition/subtraction of small number of 
particles from the MI regime. 

We have used QMC simulations to determine where 
first order phase transitions become second order. Con- 
tinuous MI-SF transitions are depicted in Fig. |3] by 
dashed mean field lines and open squares representing 
QMC results. For z = 4, the n = 1 (n = 3) MI-SF tran- 
sition becomes of second order at zJ/U = 0.20 ± 0.02 
{zJ/U = 0.133±0.02). Finally, the n = 1 Ml-to-n = 3 MI 
transition becomes second order at zJ/U = 0.040±0.008, 
setting the triple point of the phase diagram. We can 
also estimate the triple point using a three-state ap- 
proximation where we truncate the Hilbert space to the 
rii = 1,2,3 states. In this approximation the mean field 
Hamiltonian becomes: 

/-/i + 4JV^2 _4v/2j^ \ 

-4V2 Ji/j /7-2/i + 4JV^2 -AVSJiIj 
\ -AVSJtjj 3/7 + [73 - 3/i + 4 JV^2 J 



We analyze the eigenvalues of the three-state Hamilto- 
nian at /i = /ii3 and \Us\ > U. The triple point corre- 
sponds to the (J/U) TP value at which the insulator to SF 
transition becomes of second order, that is, when the cur- 
vature of E[iIj] ai ip = changes sign becoming negative. 

Solving for lim = leads to z U)^^ = 

i.e. {jj)rj.p = 0.05 for U3 = -1.5/7 and z = 4. 

First order phase transitions present in our model can 
be experimentally detected due to a loss of adiabaticity 
across the phase boundary even upon an arbitrarily slow 
ramping up or down of the optical lattice, as suggested 
in [28 or by observing the hysteretic behavior. In ad- 
dition, first order phase transitions are characterized by 
discontinuity in density profiles, a local observable easily 
accessible with state of the art techniques [29] . 

Finally we have looked at how first order phase transi- 
tions are affected by finite temperature. Strictly speak- 
ing, the MI state exists only at zero temperature. In prac- 
tice, MI features persist up to temperature T ~ 0.2/7 [30j. 
QMC results show that phase transitions are still of first 
order at temperatures of T ~ J, where MI features are 
still well defined. A more extensive study of the behavior 
of the system at finite temperature will be the subject of 
future investigations. 

Concluding, we have studied an extended version of 
the Bose-Hubbard model which includes an attractive 
3-body interaction term Us^ both at a mean field level 
and by means of quantum Monte Carlo simulations. The 
3-body term results from a universal three-body bound 
state attached to an excited threshold and can be tuned 
via an external rf field. We have found that at |/73| > /7, 
where the n = 2 lobe disappears, there exists a first or- 
der phase transition separating the n = 1 from the n = 3 
lobes which extends up to a triple point. A strong 3-body 
attraction also affects the order of the MI-SF transition. 
We have found first order transitions separating the n = 1 
and n = 3 MI from the SF. We have also analyzed the 
effect of finite temperature and found that transitions are 
still of first order at temperatures T ~ J. In the future 
it would be interesting to extensively address finite tem- 
perature properties of the system and study the case of 
repulsive three-body interaction. 
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